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Abstract 

We study neutrino oscillations and the level-crossing probability -Plsz in 
power-law potential profiles A(r) oc r n . We give local and global adia- 



O ■ baticity conditions valid for all mixing angles i? and discuss different rep- 

resentations for Plsz • For the 1/r 3 profile typical of supernova envelopes 
we compare our analytical to numerical results and to earlier approxi- 

•1— I , 

^ . mations used in the literature. We then perform a combined likelihood 

analysis of the observed SN 1987A neutrino signal and of the latest solar 
neutrino data, including the recent SNO CC measurement. We find that, 
unless all relevant supernova parameters (released binding energy, v e and 
Ua T temperatures) are near their lowest values found in simulations, the 
status of large mixing type solutions deteriorates considerably compared 
to fits using only solar data. This is sufficient to rule out the vacuum- 
type solutions for most reasonable choices of astrophysics parameters. 
The LOW solution may still be acceptable, but becomes worse than the 
SMA-MSW solution which may, in some cases, be the best combined 
solution. On the other hand the LMA-MSW solution can easily survive 
as the best overall solution, although its size is generally reduced when 
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compared to fits to the solar data only. 
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I. INTRODUCTION 



The detection of astrophysical neutrinos has played a major role in the still on- 
going process of establishing neutrino masses and mixing The propagation of these 
neutrinos from their source to the detector can be influenced by matter effects as it was 
pointed out by Wolfenstein [[| • In particular, Mikheev and Smirnov H showed that the 
flavour of neutrinos produced in the core of the sun or a supernova can be efficiently 
converted into a different one even for small mixing angles $ via "resonant" neutrino 
oscillations (MSW effect). In the limit that the density inside a star varies much 
slower than the typical distance characterizing flavour oscillations, the instantaneous 
matter states of the neutrino propagating towards the vacuum change adiabatically 
their flavour composition which can be reversed completely even for small 

The analytical study of non-adiabatic neutrino oscillations started soon after the 
discovery of the MSW effect. The leading non-adiabatic effects were calculated for 
a linear potential profile in Ref. 0] as a Landau-Stiickelberg-Zener (LSZ) crossing 
probability ||, 

^lsz = exp (-^) . (1) 
The adiabaticity parameter for a linear profile is 



7 = 
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where E is the energy, denotes the masses, and t9 the vacuum mixing angle of the 
two (active) neutrinos. Furthermore, A = 2EV = 2\f2GpN e E is the induced mass 
squared for the electron neutrino. The parameter 7 has to be evaluated at the so- 
called resonance point, i.e. the point where the mixing angle in matter is ?9 m = tt/4. 
For a linear profile, adiabaticity is maximally violated at this point. Therefore, the 
probability that a neutrino jumps from one branch of the dispersion relation to the 
other one is indeed maximal at the resonance. Since the condition $ m = tt/4 can be 
fulfilled only for a normal mass hierarchy in the case of neutrinos and for an inverted 
hierarchy in the case of anti-neutrinos, Eq. (|2|) allows the calculation of Plsz only in 
half of the parameter space of neutrino mixing. 

Later, Kuo and Pantaleone derived the LSZ crossing probability for an arbitrary 
power-law like r n potential profile. This type of profile does not only contain the 
case n « —3 typical for supernova envelopes, but also the exponential profile of the 
sun in the limit 77 — > ±00. Moreover, it allows discussing which features of neutrino 
oscillations are generic and which ones are specific for, e.g., the linear profile 77 = 1 
usually discussed. Kuo and Pantaleone found that also for arbitrary power-law like 
potential profiles the dependence on the neutrino masses and energies can be factored 
out, while the effect of a non-linear profile can be encoded into a correction function 
T 

•J m 



ftffi = exp(-^ , (3) 

which only depends on $ and n. The "adiabaticity" parameter 7 n has to be evaluated 
still at the resonance although, as we will show, it does not coincide with the point 
of maximal violation of adiabaticity (PMVA) for n 1. An unsatisfactory feature 
of Eq. (^|) is its restricted range of applicability: As in the case of a linear profile, 
the calculation of Plsz with Eq. @ requires the validity of the resonance condition 
$m = 7r/4. Several other authors have also considered non-adiabatic effects in neutrino 
oscillations, partially for arbitrary potential profiles, but their results either do not allow 
simple numerical evaluation or have also a restricted range of validity |7| . Therefore, 
it has not been possible to calculate, e.g., the survival probability of supernova anti- 
neutrinos in the quasi-adiabatic regime without solving numerically their Schrddinger 
equation. 

The purpose of this work is twofoldQ. First, we clarify the physical significance of 
the resonance point compared to the point where adiabaticity is maximally violated 
locally: we find that the crossing probability has its maximum at the PMVA and not 
at $ m = tt/4. Then we show explicitly that the product •y n J r n can be evaluated at an 
arbitrary point. We conclude that the "resonance" point has in general no particular 
physical meaning: it does not necessarily describe the point of maximal violation of 
adiabaticity nor is it necessary to calculate the adiabaticity parameter at the resonance. 
We provide a criterion that measures the global cumulative non-adiabatic effects along 
the neutrino trajectory and gives an estimate of the border between the adiabatic and 
non-adiabatic regions for a power-law profile. Moreover, we obtain an accurate formula 
for Plsz which is valid for all d and convenient for numerical evaluation. 

Second, we apply this formula to neutrino oscillations in supernova (SN) envelopes 
and in particular to the neutrino signal of SN 1987A. Performing a combined likelihood 
analysis of the observed neutrino signal of SN 1987A and the updated global set of 
solar neutrino data including the recent SNO CC measurement PJToH, we find that the 
supernova data offer additional discriminating power between the different solutions. 
We find that, unless all relevant supernova parameters (released binding energy, v e 
and v^ T temperatures) are near their lowest values found in simulations, the status 
of the large mixing solutions to the solar neutrino problem deteriorates considerably 
compared to fits using only solar data. This is sufficient to rule out the vacuum-type 
solutions for most reasonable choices of astrophysics parameters. The LOW solution 
may still be acceptable, but becomes worse than the SMA-MSW solution. In contrast 
the LMA-MSW solution can easily survive as the best overall solution, although its 
size is generally reduced when compared to fits to the solar data only. In the analysis of 



1 Some of the results of this work have been already briefly presented in Ref. ||. 



the solar neutrino data which we adopt in the present paper, the SMA-MSW solution 
is absent at 3 standard deviations if solar data only are included but may reappear once 
SN 1987A data are included and, may, in some cases, be the best combined solution. 



II. NEUTRINO EVOLUTION: RESONANCE AND ADIABATICITY 
CONDITIONS, MAXIMAL VIOLATION OF ADIABATICITY 

We consider neutrino oscillations in a two flavour scenario and label always the 
heavier neutrino mass eigenstate with "2" . Then A is positive and the vacuum mixing 
angle d is in the range [0 : tt/2]. As starting point for our discussion, we use the 



evolution equation for the medium states ip first given in Ref. [|TT 




iA m /(4£) -tf' m \ / & 



(4) 



Here, 



A m = J (A - Acos 2$) 2 + (A sin 2$) 2 (5) 



denotes the difference between the effective mass of the two neutrino states in matter, 
$ m is the mixing angle in matter with 

A sin 2$ . , 

tan2w m = (6) 

A cos 2$ -A y ' 



and $' m = d$ m /dr. Following Ref. |T|], we rewrite the evolution equation as 

(7) 
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(9) 



and 



Asin(2^ m -2^) 

sin2tf m ' 1 ' 

The traditional condition for an adiabatic evolution of a neutrino state along a 
certain trajectory is that the diagonal entries of the Hamiltonian in Eq. (|7|) are large 
with respect to the non-diagonal ones, |A m | ^> |4i?^' m |. This condition measures 
indeed how strong adiabaticity is locally violated. Therefore, the PMVA is given by 



the minimum of A m /^ m . Differentiating Eq. (||) for a power law profile, A(r) oc r n , we 
find the minimum at 



cot(2tf m - 2d) + 2 cot(2tf m ) - - [cot(2tf E 



2#) - cot(2# m )] = 0. 



For n = 1, the PMVA is at i9 m = it/ 4 for all d. Thus, in the region where the resonance 
point is well-defined, it coincides with it. In the general case, n ^ 1, the PMVA agrees 
however only for $ = with the resonance point. Finally, we recover the result of 
Ref. [12] for an exponential profile in the limit n — > ±oo. 

In Fig. [I], we show the survival probability p(r) = I^MI 2 for a neutrino produced 
at r = as z>2, together with the PMVA predicted by Eq. (|11]) and the resonance 
point for a power law profile A oc r~ 3 . The resonance condition predicts a transition in 
regions of density lower than that which characterizes the PMVA, until for $ = n/4 the 
resonance point reaches r = oo and the concept of a resonant transition breaks down 
completely. If one plots the change of the survival probability, dp(r)/dr = d\ip2{ r )\ 2 /dr, 
as function of r, cf. Fig. |], it can be clearly seen that Eq. (|TTD describes quite accurately 
the most probable position of the level crossing, while the resonance condition fails. As 
an immediate consequence, we note that in the case that the true potential profile A(r) 
is only approximately given by a power-law, its exponent should be determined by the 
region around the PMVA, not by the region around the resonance point. Finally, Fig. Q 
shows that the crossing probability becomes less and less localized near the PMVA for 
larger mixing angles 

Let us now discuss the condition for the adiabatic evolution of a neutrino state along 
a trajectory from the core of a star to the vacuum. While the condition | A m | 3> 1422$^! 
indicates whether adiabaticity is locally violated, we need now a global criterion that 
measures the cumulative non-adiabatic effects along the trajectory from $ m « 7r/2 to 

For a non-adiabatic evolution of the neutrino state we require that 
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cos^m, ip2 ~ sin^m and Eq. 
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(13) 

A/Rq, where 



We consider first the simple case of an exponential profile, |dA/dr| 
we can solve the $ m integral of the RHS of Eq. (|13| ) analytically. Then the non-adiabatic 
region obeys the following condition 

A /(#) 2 



E 



sin ■&) sin 2 (2i?) Rq 



with 
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(14) 



(15) 



This criterion agrees well with the one derived in a more intuitive way in Ref. [12 



A tan(vr/4 + fl/2) 1 

E 2sintftan$ Rq ' 1 J 

for e ps 2e'. 

In the case of a power-law like profile, A(r) = 2-EVo(r/i? ) n , the border between 
the adiabatic and non-adiabatic regions is given by 



n 

A_f f(d) 2n(2V ) l ' n y +1 

E ~ I sin 2 (2tf)(l-sintf) R~ J 

with 



[ d$ m sintf m [sin(2tf m )] 2+1/n [sin(2tf m - 2tf)] 1 ~ 1/n 

Jn/2 



(17) 
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In Fig. |H we compare the different predictions for the borders between the non-adiabatic 
and adiabatic regions for solar neutrinos with energy E = 1 MeV with the contours 
of constant survival probability (dashed lines) of the neutrino eigenstate z> 2 obtained 
by solving the Schrodinger equation (|7|). In the solar case, we can approximate the 
potential by an exponential profile with the typical solar height scale Rq = _R /1O.54. 
The difference between our general criterion (|14]) together with ( |i~5|) compared to (p~6|) 
is negligible for e' = e/2 = 1: both criteria predict very well the borderline of the 
adiabatic region. For comparison, the 90, 95, 99 and 99.73% confidence level (C.L.) 
contours for the different solutions to the solar neutrino problem from Ref. |§ are also 
shown. 

In Fig. |] we show a similar comparison for anti-neutrinos with energy E = 20 MeV 
and a profile typical for supernova envelopes, V(r) = 1.5 x 10~ 9 eV (10 9 cm/r) 3 . We 
now compare the borderlines obtained using our general criterion (|17D together with 
two different choices for /(■$): the solid line is obtained solving numerically (|T8| ) for 
n = —3, while we used (|l"5|) for the dash-dotted line. We find that already for moderate 
|n| values the function /(■#) depends rather weakly on n so that the expression ([TB]) 
obtained for n — > ±oo can be used as a reasonable approximation already for n = —3. 



III. THE CROSSING PROBABILITY IN THE WKB FORMALISM 

A. The correction functions T n 

We discuss now in detail the calculation of the leading term to the crossing proba- 
bility within the WKB formalism. The semi-classical transition probability Plsz from 
the neutrino states 1 — ► 2 under an adiabatic change of the action is mainly deter- 
mined by those x values for which Ei(x) = E 2 {x). Using the ultra- relativistic limit 
and omitting an overall phase, the WKB formula gives 
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(A-AC) 2 + (ASf 



1/2 



(19) 



where A 2 = A(C + iS) = Ae 2l1? is the branch point of A m in the upper complex x plane 
and we have also introduced the abbreviations C = cos 2d and S = sin 2d. 

We identify the physical coordinate r £ [0 : oo] with the positive part of the real x 
axis, i.e. we consider a neutrino state produced at small but positive x propagating 
to x = oo. Then, a convenient choice for A\ is to use the real part of A 2 for C > 0, 
i.e. the "resonance" point A\ = AC. However, we stress that this choice has technical 
reasons and makes sense only for C > 0: consider for instance the simplest case n — 1. 
Then both the integration path chosen and the branch cut are for C < in the half- 
plane dt(x) < 0. The physical interpretation is therefore that an anti-neutrino state 
created at small but negative x propagates to — oo. This case is however equivalent to 
a neutrino state propagating with C > in the right half-plane and therefore contains 
no new information. Thus, we expect the correction functions T n obtained with the 
integration path from AC to Ae 2l1? to be functions with period 7r/4 and to be valid 
only in the resonant region. 

We substitute first A = A(B + C), 
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Then we expand the Jacobian for a potential A(x) = Aq(x/Ro) 
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into a power series around the arbitrary point Bo. Using 

= _ ( joy /n i 
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we obtain 
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After expanding the binomial (B — Bo) m , we can solve the integral in Eq. (p0|) for 
each of the terms of (B^) separately. We obtain as our final result for general B and 
C > 0, 
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It is now clear that the choice of the point xq or Bq where the adiabaticity parameter 
7 n should be evaluated is totally arbitrary A change of B will always be compensated 
by a change in the correction function T n so that the physical observable -Plsz is 
independent of B , as it should. The choice of the resonance point B = is privileged 
only by the fact that it results in the simplest analytical expression for T n in the 
resonant region. Indeed, by requiring the invariance of Plsz against variations of Bq, 
one can write a differential equation for dJF n /di? [§]. The solution of such equation 
allows to rescale the correction functions obtained for B = to arbitrary B , 

B ) = — — F n (0, 0) . (27) 



For practical calculations however it is simpler to insert Bq = whenever possible 
(resonance case) into Eq. ( ^4|) thereby obtaining the well-known result of Ref. ||, 



Mm = 2 t f ^ 1 ) f V * ) (tan(2^)) 2m , (28) 

m=0 \ ^ m / \ 171 + 1 / 

valid only for $ < ir/8. Representing this series as a hypergeometric function, 

/n — 1 2n — 1 \ 
F n ($,0) = 2 F 1 ^__,__;2;-tan 2 (2^)j , (29) 

we can use the Euler integral representation [131 of 2 Pi, 

2 F 1 (a,b ]C] z) = T{b ^_ b) f^tt b -\l-ty- h -\l-tz)\ (30) 

as the analytical continuation for all $ G [0:7r/4]. 



B. The correction functions Q n 

In order to avoid the limited regime of validity of the formalism presented in the pre- 
vious subsection we will now present a new representation for the crossing probability 
which is valid for all We start directly from Eq. (|19|), but use now as integration path 
in the complex x plane the part of a circle of radius A centered at zero and starting from 
A\ = A and ending at A 2 = Ae 2n? namely, the end of the branch cut, see Fig. g. In 
the case of a power-law potential A = A (r/R ) n by substituting x = R (A/A o y/ n e 1,fi 
we can factor out the $ dependence of Plsz into functions Q n , 
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(32) 



(33) 



The functions Q n are well suited for numerical evaluation and always correspond to a 
neutrino state propagating in the physical part of the x plane, x > 0. Therefore, they 
have, in contrast to the T n functions, the period 7r/2 and are valid for all °d. Anti- 
neutrinos feel a potential A with the opposite sign than neutrinos. This sign change 
can be compensated by the exchange of d with 7r/2 — i?; thus the formulae derived 
for the crossing probability of neutrinos become valid for anti-neutrinos substituting 
^vr/2-0. 

Finally, we remind the reader that the Landau-Stiickelberg-Zener approach is only 
valid for small deviations from adiabaticity. Therefore, in order to cover also the non- 
adiabatic case, the WKB formula Eq. ([|) has to be replaced || by 

p _ exp (-K n g n ) - exp (-K n g' n ) 

l-expC-zc^O ' [M) 

where Q' n = Q r J sin 2 ■& for neutrinos and Q' n = Q n j cos 2 i? for anti-neutrinos, respectively. 
A similar formula holds for the T n functions, if n n Q n is replaced by 7^-" n 7r/2. 



IV. NEUTRINO OSCILLATIONS IN SUPERNOVA ENVELOPES 

A. General discussion 

The potential profile A(r) in supernova (SN) envelopes can be approximated by a 
power law with n w —3, and V(r) = 1.5 x 10~ 9 eV (10 9 cm/r) 3 |[L4|| . Since only u e 
were detected from SN 1987A and also in the case of a future galactic SN the v e flux 
will dominate the observed neutrino signal, an analytical expression for P c valid in the 
non-resonant part of the mixing space is especially useful. In the following, we will 
always consider oscillations of anti-neutrinos. 

The probability for a v e to arrive at the surface of the Earth can be written as an 
incoherent sum of probabilities, 

P- e - e = P^P* + PgPi = (1 - P c ) cos 2 d + P c sin 2 , (35) 

where Pj denotes the probability that a v e leaves the star as mass eigenstate z/j and 
Pjf the probability that Vi is detected as v e . 



In Fig. |, we compare the results of a numerical solution of the Schrodinger equa- 
tion (|7|) with the analytical calculation of P ss using the and the JF_ 3 functions. The 
latter is shown only for its range of applicability, *& > it /A. The agreement between 
the two methods using the WKB approach is again (for d > ir /A) excellent. Generally, 
these two methods agree also very well with the results of the numerical solution of the 
Schrodinger equation; there are only tiny deviations in the region A/(2E) ~ 10~ 17 eV. 

Next, we discuss the quality of the different approximations for Plsz which have so 
far been used in the literature. The most common one is to use the correct potential 
profile, A oc r~ 3 , together with T\ = 1 [15]. Since all correction functions T n go to 
1 for small mixing in the resonant region, this approximation is certainly justified for 
d — > in neutrino and for $ — > n/2 in anti- neutrino oscillations. The shape of the 
MSW triangle which is determined^] by the exponent n is correctly reproduced in this 
approximation, as the comparison of the numerical results and the linear approxima- 
tion in Fig. [7] confirms. However, this approximation becomes worse for tan 2 $ < 5 and 
breaks down for tan 2 d < 1. Therefore, it is not applicable in the phenomenologically 
most interesting region of maximal or nearly maximal mixing presently indicated by 
the solar neutrino data P, p!0] , ^6[ . Physically, the failure of the linear approximation 
15| already in the resonant region is understandable by the broadening of the cross- 



ing probability profile, cf. Fig. |2|. Since for larger mixing the crossing probability is 
determined not only by the potential near the expansion point but by a rather large, 
extended region, a linear approximation of the potential profile is not good enough. 

In order to obtain an approach well-defined for all $ another approximation was 
employed in Ref. ]T7|]. There, the expression 7 n JF n valid for an exponential profile was 
used. The resulting shape of the MSW triangle outside the region of nearly-maximal 
mixing is not correctly reproduced because of the wrong profile used, as can be seen 
from Fig. || The agreement between P c from the numerical solution of the Schrodinger 
equation and that obtained by using the exponential profile is however very good for 
tan 2 $ < 5 and a suitable choice of Aq and Rq. Therefore, this approximation is 
adequate for the discussion of anti-neutrino oscillations in the whole parameter space 
allowed by present solar neutrino experiments. Numerically, we have found the best 
agreement between the two approaches for R^ p w 100 x i?Q° wer , where we have fixed 
A^ p = AQ° wer . This normalization of the scale factor R ensures that both potentials, 
A = y4 exp(— r/i?Q Xp ) and A = A (r / R^ oweT )~ 3 , give rise to the same scale height, 
[d In A(r) / dr]" 1 , in the region around the PMVA for the values of $ of interest. Figure § 
shows very good agreement between this approximation and the numerical results 



2 Note that, according to Eq. (Q), the hypotenuse of the MSW triangle should have the slope 
n/(n + 1) in a ln(A)-ln(tan 2 plot. 



except for tan 2 $ > 3 and A/E > lO" 10 eV 2 /MeV. 

Reference |J gave an additional approximation for JF_ 3 which consists of an expan- 
sion of the exponential profile in the parameter l/n, 

F n (d) = (1 - tan 2 ??) {1 - — [ln(l - tan 2 tf) + 1 



1 + tan 2 ■&) 



ln(l + tan 2 $)] + ...}. (36) 



tan 2 ■& 

This approximation agrees indeed very well already for n = —3 with our numerical 
results in the resonant region, cf. Fig. [7J Only near the limit of its range of applicability, 
tan 2, # = 1, sizeable deviations can be seen. 

Finally, we have also checked how deviations of the true SN progenitor profile V(r) 
from an 1/r 3 profile may affect our results for Pg g . Realistic progenitor profiles differ in 
two aspects from the simple 1/r 3 profile assumed. First, the outer part of the envelope 
has an onion like structure, and its chemical composition, Y e (r), and thus also V(r) 
changes rather sharply at the boundaries of the various shells. Second, the density drops 
faster in the outermost part of the envelope, becoming closer to an exponential decrease. 
We have calculated numerically P gg using profiles for three different progenitor masses 
(11, 20 and 30 M ) from Woosely [18] and one by Nomoto ][9] for the progenitor of 
SN 1987A. We have found that P gg is well approximated in the non-resonant part by 
our analytical results for the 1/r 3 profile, independently of the details of the progenitor 
profile. In contrast, P gg depends strongly on the details of the progenitor profile in 
the resonant region. In particular, there are rather strong deviations between a 1/r 3 
profile and the profiles of Refs. [18,191. As an example, we compare in Fig. || the 
P S e calculated numerically for the 2OM profile of Ref. |18 with the analytical results 
for our standard SN profile V{r) = 1.5 x 10~ 9 eV (10 9 cm/r) 3 . We conclude that a 
numerical solution of the Schrodinger equation (|7]) should be performed in the resonant 
region, using a realistic profile for the particular progenitor star considered. However, 
a 1/r 3 profile together with the LSZ crossing probability is sufficient for the analysis of 
antz-neutrino oscillations in the phenomeno logically most interesting region tan 2 $ < 5. 
Note that it is this region where the other methods based on the splitting of Plsz into 
the adiabaticity factor 7 plus the functions T fail. 



B. Likelihood analysis of the SN 1987A neutrino signal 



In Ref. [17] we performed a likelihood analysis PUj of the neutrino signal of 
SN 1987A observed by the Kamiokande II and the 1MB experiments |2l] using the 
expression ^i n T n valid for an exponential potential profile for Plsz- in order to describe 
the SN 1987A neutrino signal we need the evolution of a typical SN anti-neutrino in the 
region of parameters presently indicated by the solar neutrino problem ||9|,[T0|] . Most of 



the presently allowed oscillation solutions of the solar neutrino problem (LMA-MSW, 
LOW-QVAC, VAC, "just-so") require large neutrino mixing. The small mixing an- 
gle solution, called SMA-MSW, appears in most analyses only at about 99.9% C.L. 
The precise value slightly depends on arbitrary details of the statistical analysis, as 
can be seen by comparing different global analyses of the solar data. So far, only 
the SuperKamiokande collaboration can perform a full analysis of its signal and back- 



ground ||22|| . However, the published data allow to approximate reasonably well the full 
SuperKamiokande results so as to be used in global analyses. We have here adopted 
the values of the X©C$j A) grid corresponding to Fig. 6b of Ref. ||. 

The description of the anti-neutrino conversion is trivial both in the cases of com- 
plete adiabaticity or full non-adiabaticity. In the first case there is, of course no-level 
crossing, while in the second the standard vacuum treatment applies. In either case 
the exact form of the potential profile used has no effect. As we have seen in Figs. § 
and |] full adiabaticity or non-adiabaticity occur for an exponential and a power-law 
like profile with n = —3 for all solutions to the solar neutrino problem, except for the 
case of LOW-QVAC solution. Therefore for the former solutions the main part of the 



results and conclusions of Ref. [17] remains unchanged. However, for the LOW-QVAC 
a new analysis is in order. 

Here we have repeated therefore the likelihood fit of the neutrino signal from 
SN 1987A in the tan 2, #-A plane for a more realistic 1/r 3 profile together with the 
expression for the crossing probability but used otherwise exactly the same ap- 
proach as in Ref. |17]]. In particular, we work in the framework of oscillations between 
two active neutrino flavours. Our 2-flavour approximation is valid for the case of nor- 
mal mass hierarchy m 3 ^> m x ~ m 2 in view of the Chooz limit [23] on U13. In the 
alternative inverse mass hierarchy case mi ~ mi 3> m 3 , the third neutrino do not 
significantly affect supernova oscillations if sf 3 < few x because then the "heavy" 
resonance is non-adiabatic. With the inverted spectrum, larger values of S13 are already 



strongly disfavoured by SN 1987A pi |. 

The ratio of the likelihood functions L for different hypotheses measures the degree 
to which the experimental data favour one hypothesis over the other. In order to decide 
how strong the large mixing solutions are disfavoured with respect to the no-oscillation 
case (or to SMA-MSW oscillations, that negligibly affect supernova v e ) for a certain 
choice of astrophysical parameters {a}, we consider the ratio 

Ra(fi,*)= r ai ^ A) ( 37 ) 

*-a,NO-OSC 

as function of $ and A. The likelihood function 

£ a (tf, A) oc exp \ - J n a (E) dEj JJ n a (E h tf, A) (38) 

tests the hypothesis that a prescribed neutrino fluence F Pe leads to the observed exper- 
imental data Ei with probability distribution n a (E, A). The connection between the 



observed data set Ei and the emitted neutrino fluence F has been described already in 



detail in the literature, cf. Ref. 17 



The astrophysical parameters a used are the released binding energy E^, the average 
energy {E Dc ) of v e neutrinos, and the average energy (E Ph ) = r(E 9e ) of P MT neutrinos. 



For a more details see Ref. |17| . We consider mainly two sets of values for {E^, (E Pe )}: 
one corresponds to the lowest values found in simulations, E^ = 1.5 x 10 53 erg and 
(E Dc ) = 12 MeV, and the second one is with E h = 3 x 10 53 erg and (E Pa ) = 14 MeV 



nearer to average values found [25 



In Figs. 10 and 11 we show the likelihood ratio ln(i?) as function of tan 2 $ and A 
relative to the NO-OSC hypothesis. The contours of constant likelihood shown cor- 
respond to ln(i?) = — 1, —2, —3, —5, —10, —15, —20, unless otherwise indicated. There 
are two differences compared to the same plots of Ref. First, since in Ref. |TTJ 

the scale height Rq of the exponential profile was not optimized to reproduce best 
a power-law profile with n = —3, matter effects in the SN envelope became impor- 
tant at 10~ 7 -10~ 8 eV 2 while here we find that they already start being important for 
few x 10~ 9 eV 2 . Second, the slope of the MSW triangle in the dark side tan 2 $ > 1 has 
now the correct value. The changes in the regions favoured by the solar neutrino data 
are however marginal. 



C. Combined global analysis of solar and SN 1987A neutrino data 

In Figs. [H] and [Ll|, we have simply superimposed the CL. contours obtained by 
analyzing solar and SN data separately. Since the two data sets are statistically inde- 
pendent and functions of the same two fit parameters, they can be trivially combined, 

X 2 OT (tf, A) = A) + xU$, A) . (39) 

We can then obtain new CL. contours which are defined relative to the minimum of 
Xtot- m ^e following, we will use the x©($> A) grid calculated in Ref. || for the 
solar data and xIn = ~ 2 A ) with £(0, A) as defined in Ref. |7| for the SN 1987A 



data. We consider the astrophysical parameters as known and minimize only the two 
parameters d and A. Hence, the confidence levels are always calculated with respect 
to two fit parameters. 

In Figs. p!2|-[T6| we show the 90, 95, 99 and 99.73% CL. contours of the combined 
fit of solar and SN 1987A data for different astrophysical parameters. We also show 
the best-fit point of the solar neutrino data alone (dot) and the best-fit point of the 
combined data set (star). The values chosen in Fig. [12] correspond to a rather represen- 
tative set of astrophysical parameters, namely E^ = 3 x 10 53 erg and (Ep c ) = 14 MeV. 
In this case, the impact of the SN 1987A data on the standard solutions to the solar 
neutrino problem is rather dramatic: the LOW-QVAC and VAC solutions disappear 



for all assumed r values. We find that they are excluded at 99.98% and 99.99% respec- 
tively, even for r = 1.4. Moreover the size of the LMA-MSW solution decreases, with 
increasing r. The part of the LMA-MSW solution which is most stable against the ad- 
dition of the supernova data corresponds to the lowest A and tan 2 $ values, since these 
are favoured by Earth matter regeneration effects. On the other hand the SMA-MSW 
region re-appears extending, for increasing r, as a funnel towards the VAC solution 
along the hypotenuse of the solar MSW triangle. The combined best-fit point moves 
from the LMA-MSW region for r = 1.4 to the SMA-MSW solution for r = 1.7 and 2. 
Comparing the C.L. contours for the different solutions to the solar neutrino problem 
with the border between the adiabatic and non-adiabatic region in Fig. we find that 
the SMA-MSW solution lies for anti-neutrinos always in the adiabatic region, while it 
is in the non-adiabatic region for neutrinos (dotted line). 

In Fig. |13, we have chosen astrophysical parameters corresponding to the lower 
limit of the range found in simulations, = 1.5 x 10 53 erg and (E Pc ) = 12 MeV. 
Even in this case, the vacuum type solutions to the solar neutrino problem are severely 
disfavored, and exist very marginally only at the 3a level for r — 1.4. In contrast, the 
LOW solution can still exist at 90% level for the most favorable choice of astrophysical 
parameters. Finally the LMA solution remains for all of the above three r choices 
and, most importantly, for such choices of astrophysical parameters the best combined 
global fit-point remains within the LMA solution. 

In order to analyse the stability of the LMA solution against larger values of the 



astrophysical parameters, we have considered in Fig. the following case, = 
3 x 10 53 erg, and (E Pc ) = 16 MeV. It can be seen how even for this choice the LMA- 
MSW solution is allowed at 99% C.L unless (E Ph ) becomes uncomfortably large (r = 2), 
in which case no large mixing solutions remains at 3a. 



To elucidate more the importance of the different parameters, we show in Figs. [L5 



and the confidence contours for low E^ combined with an average value of (E Pc ) 
and vice versa. Comparing the two set of figures, one recognizes that a low value of E^ 
can be more important for the "goodness-of-fit" of the LMA-MSW region than a low 
value of (E Pc ). In the case of small energy release, E^ = 1 x 10 53 erg, the combined 
best-fit point remains even for r = 2 and (E Pe ) = 14 MeV in the LMA-MSW region. 

In Fig. |T7] we illustrate in a global way the relative status of various oscillation solu- 
tions to the solar neutrino problem and the qualitative impact of adding the SN 1987A 
data. For each A value we have optimized the x 2 with respect to The solid curve 
indicates the \ 2 of the various oscillation solutions to the solar neutrino problem. The 
non-solid curves correspond to the case where the SN 1987A data are included. The 
dash-dotted line is for E h = 3 x 10 53 erg, r = 1.4 and (E Pc ) = 14 MeV. The dashed 
line is for E h = 3 x 10 53 erg, r = 1.4 and (E Pc ) = 12 MeV. The dotted line is for 
£b = 3 x 10 53 erg, r = 1.7 and (E Pc ) = 14 MeV. Here we have adjusted an arbitrary 
constant which appears when combining the minimum likelihood-type SN 1987A anal- 



ysis with the solar x 2 data analysis in such a way that the SMA solution gets unaffected 
by the SN 1987A data. One notices that the effect of adding SN 1987A data is always 
to worsen the status of presently preferred oscillation-type solutions. Within each such 
curve one can compare the relative goodness of various solutions, however different 



curves should not be qualitatively compared. In Fig. [18| we repeat this procedure after 
optimizing over A and displaying the result with respect to the mixing parameter 

In contrast to earlier investigations, we conclude that the observed neutrino signal 
of SN 1987A is compatible with the LMA-MSW solution, especially for "standard" 
values of the neutrino energies together with small values of E\> Smaller apparent 
values for are e.g. expected if SN 1987A was not spherically symmetric and the axis 
of maximal emission of neutrinos was not pointing to the Earth. We note also that, 
although the size of the LMA-MSW region in the combined fit may be substantially 
reduced, the position of its local best-fit point is rather stable: the inclusion of the 
SN 1987A data drives the local best-fit point to only slightly smaller values of d and 
A. While there is no significant conflict between the solar and SN 1987A data for the 
LMA-MSW solution, the case of the other large mixing solutions is different: even the 
LOW solution which always gives a better fit than the vacuum-type solutions requires 
a conspiracy of all astrophysical parameters to fit well the combined data. 



V. SUMMARY AND DISCUSSION 

We have discussed non-adiabatic neutrino oscillations in general power-law poten- 
tials A oc x 11 . We found that the resonance point coincides only for a linear profile with 
the point of maximal violation of adiabaticity. We presented the correct boundary 
between the adiabatic and non-adiabatic regime for all d and n. As a new method 
to calculate the crossing probability also in the non-resonant regime we proposed a 
different splitting of In P c into a $ independent part n n and new correction functions 
Q n which have a simple integral representation for all •& G [0:7r/2]. 

As an important application for supernova neutrinos, we considered the case n = 
—3 in detail. Comparing the different approximations used in the literature with 
exact numerical results, we found that all of them fail in some part of the tan 2 
A plane. While the use of T\ = 1 together with A oc r -3 describes correctly the 
crossing probability for small mixing in the resonant region, the errors become larger 
for larger mixing until this approximation fails completely in the non-resonant region. 
The correction function JF^ used for n = — 3 describes quite accurately the most 
interesting region of large mixing as well as the non-resonant region, but does not 



3 A deviation of equipartition, i.e. a smaller luminosity of v^ T than of u e , has a similar effect. 



reproduce the correct shape of the MSW triangle. In contrast, the results using the 
JF_ 3 and £7_3 function agree very well with the numerical results for a 1/r 3 profile in 
the resonant region and for all respectively. 

Performing a combined likelihood analysis of the observed neutrino signal of 
SN 1987A and solar neutrino data including the recent SNO CC measurement, we 
have found that the supernova data offer additional discrimination power between the 
different solutions of the solar neutrino puzzle. Unless all relevant supernova param- 
eters lie close to their extreme values found in simulations, the status of the large 
mixing solutions deteriorates, although the LMA solution may still survive as the best 
fit of solar and SN 1987A data for acceptable choices of astrophysical parameters. In 
particular, SN 1987A data generally favour its part with smaller values of d and A. 
In contrast the vacuum or "just-so" solution is excluded and the LOW solution is 
significantly disfavoured for most reasonable choices of astrophysics parameters. The 
SMA-MSW solution is absent at about the 3<r-level if solar data only are included but 
may reappear once SN 1987A data are added, due to the worsening of the large mixing 
type solutions. 

Finally, we speculate about the possible impact of the SN 1987A data on the solar 
neutrino problem in the near future. If the solution of the solar neutrino puzzle lies 



in the LMA-MSW region, the KamLAND experiment |26[ will determine *& and A 
accurately [[27] (unless A > 2 10 _4 eV 2 ). However KamLAND cannot distinguish t? 
from 7r/2 — -d. Both solar and SN 1987A data can make this distinction, and both 
favour d < n/A: in the solar case because it gives a MSW resonance in v e oscillations; 
in the SN 1987A case because it gives no MSW resonance in v e oscillations. On the other 
hand, the knowledge of the neutrino mixing parameters would change the emphasis of 
studies of the SN 1987A neutrino signal from particle physics to astrophysics and help 
to shed new light on the allowed range of the astrophysical parameters. 

Last but not least, let us mention that one should not forget a fundamental differ- 
ence between the SN 1987A and the solar fits. In the solar case, a well-tested standard 
solar model exists, whose errors in the various quantities from astro- and nuclear physics 
can be estimated and are accounted for in the fit. In contrast there is no "standard 
model" for type II supernovae and therefore also no well-established average values and 
error estimates for the relevant astrophysical parameters. Nevertheless, experimental 
data about SN 1987A exist, and it is worth studying their consequences. 
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FIG. 1. Survival probability p(r) = \tp2(r)\ 2 as function of r for a neutrino produced 
at r = as z/j. The point of maximal violation of adiabaticity predicted by Eq. (|ll|) is 
indicated by a dot, while the resonance point for a power law profile A = 2EV (r/R )- 3 with 
A/E = 2 x 1CT 10 eV 2 /MeV, V = 1.5 x 1CT 9 eV and R = 10 9 cm is denoted by a star. 
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FIG. 2. Change of the survival probability dp(r)/dr of a neutrino produced at r = as 
v% together with the point of maximal violation of adiabaticity (dot) and the resonance point 
(star) for a power law profile A oc r~ 3 . The height of the different curves is rescaled. 
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FIG. 3. Contours of constant neutrino survival probability (dashed) together with the 



borderlines Eq. (14) (solid) and (|^) (dash-dotted) between adiabatic and non-adiabatic 
regions for the solar density profile. The 90, 95, 99 and 99.73% C.L. contours for the different 
solutions to the solar neutrino problem are also shown. 
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FIG. 4. Contours of constant anti- neutrino survival probability (dashed) together with 



the borderlines Eq. (|i7|) between adiabatic and non-adiabatic regions using Eq. ( |18|) (solid) 
and ( [To]) (dash-dotted) for /(#) for the SN profile given in the text; the dotted line shows the 
borderline for neutrinos. The 90, 95, 99 and 99.73% C.L. contours for the different solutions 
to the solar neutrino problem are also shown. 
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FIG. 5. Analytical structure of the function A m together with the different integration 



paths used in Sec. MB for the special case n = 1. The branch cut is shown by a wavy line, 
poles by crosses; an additional pole appears for n < at x = 0. 
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FIG. 7. Comparison of the contours of constant survival probability P ss calculated 
numerically (solid lines) and with the linear approximation J~\ = 1 (dashed lines) and the 
approximation Eq. ([?]) (dotted lines), respectively; all for A oc r -3 as given in the text. 
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FIG. 8. Comparison of the contours of constant survival probability Pgg calculated nu- 
merically (solid lines) for A oc r~ 3 and with the exponential approximation A cx exp(— r/i?Q Xp ) 



(dashed lines) with = 100 x ,RP ower . 
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FIG. 9. Comparsion of the contours of constant survival probability P ss calculated 
numerically for a M = 2OM progenitor star from Ref. [17] (dotted lines) and calculated for 
A oc r -3 with the Plsz approximation (solid lines). 
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FIG. 10. Likelihood ratio ln(i?) relative to the NO-OSC hypothesis, as function of tan 2 d 
and A/eV 2 for r = (E Ph )/(E Pc ) = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom) together 
with the different solutions to the solar neutrino problem. All figures for = 1.5 x 10 53 erg 
and (E Pc ) = 12 MeV. 
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FIG. 11. Likelihood ratio ln(i?) relative to the NO-OSC hypothesis, as function of tan 2 d 
and A/eV 2 for r = (E Ph )/(E Pc ) = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom) together 
with the different solutions to the solar neutrino problem. All figures for E^ = 3 x 10 53 erg 
and (E Pc ) = 14 MeV. 
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FIG. 12. The 90, 95, 99 and 99.73% C.L. contours of the combined fit of solar and 
SN 1987A data (coloured/grey) together with the contours of the solar data alone (solid 
lines); for r = (Ep h ) / (Ep c ) = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom). All figures for 
B b = 3 x 10 53 erg and (E Pc ) = 14 MeV. 
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FIG. 13. The 90, 95, 99 and 99.73% C.L. contours of the combined fit of solar and 
SN 1987A data (coloured/grey) together with the contours of the solar data alone (solid 
lines); for r = (Ep h ) / (Ep c ) = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom). All figures for 
E h = 1.5 x 10 53 erg and (E Pc ) = 12 MeV. 
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FIG. 14. The 90, 95, 99 and 99.73% C.L. contours of the combined fit of solar and 
SN 1987A data (coloured/grey) together with the contours of the solar data alone (solid 
lines); for r = (Ep h ) / (Ep c ) = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom). All figures for 
B b = 3 x 10 53 erg and (E Pc ) = 16 MeV. 
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FIG. 15. The 90, 95, 99 and 99.73% C.L. contours of the combined fit of solar and 
SN 1987A data (coloured/grey) together with the contours of the solar data alone (solid 
lines); for r = (Ep h ) / (Ep c ) = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom). All figures for 
E h = lx 10 53 erg and (E Pc ) = 14 MeV. 
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FIG. 16. The 90, 95, 99 and 99.73% C.L. contours of the combined fit of solar and 
SN 1987A data (coloured/grey) together with the contours of the solar data alone (solid 
lines); for r = (Ep h ) / (Ep c ) = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom). All figures for 
B b = 3 x 10 53 erg and (E Pc ) = 12 MeV. 
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FIG. 17. The solid curve indicates the x 2 of the various oscillation solutions to the solar 
neutrino problem. The non-solid curves illustrate the effect of adding the SN 1987A data, 
which worsens the status of large mixing-type solutions. See text for explanation. 




FIG. 18. Same as in Fig. 17 but displayed with respect to see text for explanation. 



ADDENDUM: COMBINED ANALYSIS OF SN 1987A AND SOLAR DATA 
AFTER SNO'S NC AND DAY/NIGHT RESULTS 



In this addendum, we update our analysis including the new NC and day /night 
measurements of SNO [TJ as well as the latest Super-K ]2[], SAGE || and GNO re- 
sults . As a result of the precise determination of the 8 B flux by SNO together with 
the absence of spectral distortions in Super-K, the solar data disfavour now SMA-MSW, 
VAC, and the LOW solutions more strongly than before. To show the impact of the 
new measurements, we have performed a new combined \ 2 analysis of the SN 1987A 
and the solar data, using the same method as described in section [IV Cj . The x© were 
taken from Ref. ||, where also more details of this analysis can be found. 

We consider again the astrophysical parameters (E^, (E Pc ), r and implicitly the 
assumption of equipartition) as given a priori, and use only $ and A as fit parameters. 
Since present simulations with improved treatment of neutrino transport and micro- 
physics (inclusion of energy transfer by recoils, nucleon bremsstrahlung, and flavour- 
changing neutrino annihilations v e v e —>■ ^ iT z^ iT ) find much smaller differences between 
the spectra of v e and z/ M)T , we use now only the lowest value for the ratio r of the z/ M r 
and v e temperatures considered previously, r = 1.4 ||. There are two main conse- 
quences of low r values: First, the weight of the SN 1987A data in the combined fit is 
reduced. Second, deviations from equipartition which are naturally expected become 
more important and should not be neglected in the future. 

In Fig. we show the 90, 95, 99 and 99.73% C.L. contours of the combined fit 
for different v e average energies, (E Pc ) = 12, 14 and 16 MeV. In all cases, we use 
i?b — 3 x 10 53 erg and r = 1.4; also, we keep the assumption of equipartition to allow a 
simple comparison with the previous figures. The best-fit point of the combined data 
set lies now for all values of (E Pc ) in the LMA-MSW region, very close to the solar 
best-fit point. Both the new solar data and the small value of r strengthen the status of 
LMA-MSW as the leading solution to the solar neutrino problem. Only the regions of 
the LMA-MSW solution with large A or tan i? > 1 — where Earth regeneration effects 
are less important — are inconsistent with the combined data set. As in our previous 
analysis, we find that the other large mixing solutions, LOW, quasi- VAC and VAC, 
are disfavoured by the SN 1987A data. The LOW solution appears for small values 
of (E Pe ), but only at the 99.73% confidence level. On the other hand, the SMA-MSW 
solution appears at 99% C.L. if one goes to large values of (E 9e ). 

As a concise way to illustrate these results, we show in Fig. |20| the \ 2 in terms of 
only one variable, A or tan 2 $, respectively. The dash-dotted line is for r = 1.4 and 
(E Pe ) = 14 MeV, the dashed line for r = 1.4 and (E Pe ) = 12 MeV, while the dotted line 
shows x 2 f° r t = 1-7 and (E Pc ) = 14 MeV. Finally, note that the arbitrary constant 
which appears when the maximum-likelihood analysis of SN 1987A is combined with 



the solar x 2 analysis has again been adjusted in such a way that no oscillations (and 
consequently the SMA-MSW solution) are unaffected by the SN 1987A data. 

In summary, we find that the combined effect of current SN simulations favoring 
low values of r, together with the large x% difference between the LMA-MSW solution 
and its competitors, result in a very comfortable position of the LMA-MSW solution 
as the leading solution to the solar neutrino problem. 
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FIG. 19. The 90, 95, 99 and 99.73% C.L. contours of the combined fit of solar and 
SN 1987A data (coloured/grey) together with the contours of the solar data alone (solid 
lines); for (E Pc ) = 12 MeV (top), (E Pc ) = 14 MeV (middle) and (E Ve ) = 16 MeV (bottom). 
All figures are for E b = 3 x 10 53 erg and r = (Ep h )/(Ep a ) = 1.4. 
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FIG. 20. The solid curve indicates the x 2 of the various oscillation solutions to the solar 
neutrino problem as a function of A (top) and tan 2 i? (bottom) . The non-solid curves illus- 
trate the effect of adding the SN 1987A data, which worsens the status of large mixing-type 
solutions. See text for explanation. 





